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I.  INTRODUCTION 


The  inflationary  universe  model,  in  which  expansion  Is  accelerated  (R  >  0), 
is  expected  to  occur  in  the  early  universe  containing  matter  in  the  form  of 
bare  quantum  fields,  the  expansion  being  exponential  (R  exp[kt])  [1].  Power- 
lav  inflation,  in  which  the  expansion  scale  factor  obeys  a  power-law  relation 
with  the  time,  is  also  a  possible  result  of  those  physical  processes  in  the 
early  universe.  It  is  Important  because  it  too  will  solve  such  cosmological 
problems  as  horizons,  homogeneity,  and  flatness  [2]. 

Further,  Inflation  is  important  because  it  is  thought  that  it  can  solve 
the  problem  of  the  apparent  large-scale  isotropy  of  the  universe  [3].  This  is 
so  because  its  presence  minimizes  the  acceleration  produced  by  the  existence 
of  a  cosmological  constant,  which  lsotroplzes  the  universe  [4].  There  is, 
however,  the  model  in  the  Einstein-Cartan  gravitational  theory  with  the  Ray- 
Smalley  improved  stress-energy-momentum  tensor  (SEMT)  with  spin.  Gasperini 
has  proven  that  for  the  case  of  an  rms  spin  density  [5].  However,  it  will  be 
shown  that  the  inflationary  epoch  occurs  because  of  the  density  of  spin  angu¬ 
lar  kinetic  energy,  which  is  a  local  quantity  dependent  on  the  spin,  and 
therefore,  it  is  not  necessary  to  resort  to  an  rms  expectation  value  of  the 
spin  density  operator  to  generate  the  spin  terms  necessary  to  induce  infla¬ 
tion.  The  model  used  is  a  sample  anisotropic  (Blanchi  Type  I)  cosmological 
model  with  shear,  but  with  vanishing  spatial  curvature  (Euclidean  model). 

Formal  solutions  of  the  Einstein-Cartan  equations  and  of  the  fluid  equations 
of  motion  are  exhibited  and  shown  to  lead  to  conditions  producing  an  infla¬ 
tionary  epoch  in  the  very  early  universe.  The  Inflation  is  due  to  the  angular 
kinetic  energy  density  of  the  spin.  The  shear  is  not  effective  in  preventing 
or  damping  the  inflation  in  the  model.  The  model  further  brings  out  the 
relationship  between  the  fluid  vorticity,  the  angular  velocity  of  observers' 
lnertlally-dragged  and  Fermi -transported  reference  tetrads,  the  precession  of 
the  principal  axes  of  the  ellipsoid  of  the  shear  rate,  and  the  torsion.  It 
concludes  with  comments  of  further  work  and  suggestions  for  new  investigations 
of  this  and  related  models. 

In  Section  II,  the  basic  equations  of  the  model,  following  the  pattern  of 
Blanchi  I  spacetime,  are  given,  and  its  behavior  is  shown  in  a  Rlemann-Cartan 
spacetime  in  Section  III.  The  conclusions  are  presented  in  Section  IV. 


II.  FORM  OF  THE  MODEL 


In  this  paper,  che  metric  chosen  has  the  following  form,  as  used 
by  Mlsner  [6] ; 

ds*  *  -df  +  e-a  e*^  dx*  dx^  (1) 

where  a  is  a  scalar  function  of  time  and  exp[B]  a  traceless,  3x3  matrix, 
also  a  function  of  time.  Following  the  method  of  differential  forms,  Cartan's 
first  equations  is  written,  connecting  the  basis  forms  with  Torsion  Syv  a,  as 

.at,  a  _  v  au»  ^  .-x 

du  +  ui  A  u  *  jS  u)  A  u)  (2) 

v  uv 

where  the  torsion  is  defined  as 


S 


(3) 


Here,  the  torsion  is  defined  as  the  true  antisymmetric  portion  of  the  affine 
connection  (non-zero  in  a  holonomic  frame),  with  the  C's  the  antisymmetric 
portion  (if  any)  due  to  the  choice  of  tetrads.  Choosing  a  basis  one-form  set 
which  diagonalizes  the  above  metric  and  puts  it  into  Minkowski  form,  gives: 


i 

iii 


U1 


dt 


a  3  ,  j 

e  e  ,  .  dxJ 


(4a) 

(4b) 


In  a  Blanchi  Type  I  cosmology,  all  the  spatial  structure  constants  are 
zero.  Thus, 


jk 


0 


(5) 


The  tetrads  have  the  following  properties:  The  capital  Latin  indices  are  used 
to  refer  to  anholonomic  coordinates,  and  Greek  Indices  to  holonomic  coordi¬ 
nates.  Thus, 


EAu  V 


*AB 


EA  E,  »  n 
u  Av  uv 


(6a) 

(6b) 


Following  Ray  and 
spin-density  suu, 
is  written  as 


Smalley  [7],  an  expression  giving  the  orientation  of  the 
angular  momentum  wyv  of  the  tetrads,  and  the  Improved  SEMT 


s  ■  k(x)  (E1  E2  -  E1  E2  ). 
uv  U  v  V  u 

”uv  •  *!»<**„>  EAu  -  D(£*„)  EAu] 


(7a) 


raS 

Rav-Smalle', 


0(l  +  £  +  P/0) 


+  P 


a3 


(7b) 


Using  tetrads  consistent  with  the  choice  of  Bianchi  Type  I  (Ref.  8,  page  110) 
structure  in  Eq.  (7a),  the  non-zero  components  of  syV  are 


Again,  following  Ray  and  Smalley  [7],  the  trace-free  (proper)  torsion  saBU* 
which  is  defined  as 

S  a  =  S  a  +  (2/3)6?  S  .  3  (9' 

MV  pv  [p  v]f}  ^  ' 

is  related  to  the  spin  density  syv  by  the  relationship 

s  a  =  i<ps  ua  .  Cio; 

MV  *  H  pv 

where  K  *  8irG,  and  G  is  the  gravitational  constant.  For  purposes  of  this 
paper,  a  co-moving  frame  (ua  s  6a)  with  normalized  four-velocity  (uaua  =  -1) 
is  used.  Thus,  the  non-zero  components  of  trace-free  torsion  are: 

S12°  “  i<ps12u°  a  _S21° 

which  gives  a  relationship  between  the  torsion  and  proper  torsion: 


12  12 


itcpk(x)u°  =  -S. 


These  values  of  torsion  will  be  used  in  the  model  under  consideration. 


Proceeding  with  the  calculation,  the  connection  two-forms  are  calculated 
using  Cartan's  first  equation.  In  this  paper,  the  following  notation  for 
derivatives  are  used:  The  overdot  indicates  partial  differentiation  with 
respect  to  time,  as  in 


3/3c  A. .  or  (A. .  ) 

ij  .  .  .  ij  .  . . 


(13a) 


while  the  directional  covariant  derivative  along  the  four  veloctiy  is  indi¬ 
cated  by 


D(A. .  )  -  A. .  u 

ij . . .  ij . . . ;a 


(13b) 


The  quantity 


which  occurs  when  doing  the  computations  involved  in  Cartan's  First  Equation, 
is  split  into  a  symmetric  and  antisymmetric  part,  the  former  related  to  the 
shear,  the  latter  related  to  the  twist  of  the  congruence  of  the  normals  to  the 
homogeneous  hypersurfaces.  These  are  written 

-  e3nUI  *  (15a) 


(15b) 


The  connection  two  forms  are  summized  in  Table  1  below. 


TABLE  1.  Connection  Two-Forms 


(i  sjk  +  5jk  +  lsjk0)  "k 


5jk  +  V  “ 


-?jk  dt 


u)°  0 

o 


Using  the  connection  two-forms  summarized  above,  Cartan's  Second 
Equation  Is  computed  as 

■JO  .*  a,  -a,  *  v 

°3  dw3  +  A  *3 

-  iR%  UU  A  ^  ,u 

3uv  (1( 

Doing  the  computations  Involved  In  Eq .  (16)  are  somewhat  tedious,  and 
only  the  results  are  stated.  The  non-zero  components  of  the  Rieraann  tensor 
are  (k*j  -  o£j  +  t*j); 


a  +  o..  +  (a  +  5).,  +  [5  •  ?1,. 


(17a) 


R1  »  R1  -  R°  -  0 
ooo  ojk  ooo 


(17b) 


R° .  .  -  -R°. . 

10J  1JO 


+*u  +  <;  +  5)u  *  15  •  ;iij 


+  i(S.  .°  +  oS. ,°  +  £, .S. lo 

ij  lj  i 


x..  S.t0)  (17c) 

ll  j 


4 


(17d) 


^Jkl  ’  (i  6\  +  5V  (i  5jt  +  V  '  U  4\  +  5V  <“  5jk  +  5jk> 


+  (i  «\  +  5V  -  sjk°  ^  5\  +  51,)] 


Using  the  Reimann  tensor  components,  shown  above,  the  Ricci  tensor  com¬ 
ponents  RuU  are  then  calculated  by  contracting  the  first  and  third  indices 
of  the  Riemann  tensor.  The  non-zero  components  of  the  Ricci  tensor  are 


01 


R.  =  0, 
10 


R° 


.  7 


3  2  +  3  (a)“  +  a . .  , 

ij 


(18a) 


R-.  =  (a  3°.  +  a..)  »  3  a  (a  3..  +  j.  .) 
O  ij  tj  ij  ij 


[5 


T-] 


+  its.  0  +  3is..°  -  2r.  r .  S  . , 1 

ij  ij  k[i  j] 


(18b) 


The  curvature  scalar  R  is 


R  *  6  a  +  12  a  +5  jij  (19) 

ij 

The  16  components  of  the  Einstein  tensor  are  defined  by 

G  -  R  *  R  ,  (20) 

UV  pv  °pv 

and  summarized  in  Table  2.  For  brevity  sake,  the  additional  terms  have  been 
written  explicitly  only  due  to  torsion.  The  Einstein  tensor  which  is  obtained 
in  the  regular  GR  theory  is  listed  only  symbolically  (see  Misner  [6]  for  regu¬ 
lar  GR  components). 
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TABLE  2.  Summarization  of  the  Einstein  Tensor  Components 


u 

V 

v 

(Unprimed  G's  are  standard 
results  from  GR) 

0 

0 

Goo 

0 

i 

0 

(i  =  1,2,3) 

i 

0 

0 

(i  »  1,2,3) 

1 

1 

G11 

1 

2 

G12 

+  iisl2°  +  3 

i  su»] 

1 

3 

G13 

+  ^T23  S12 

2 

1 

G21 

+  J[S21°  +  3 

i  S21°] 

2 

2 

G22 

2 

3 

G23 

+  »13  S21° 

3 

1 

G31 

*T23  S12 

3 

2 

°32 

iT13  S21 

3 

3 

G33 

To  calculate  the  Ray-Smalley  tensor  components  in  this  model,  the  affine 
connections  will  be  needed,  which  are  contained  implicitly  in  Table  1.  Using 
the  relationship  that  waY  -  raSy  determines: 


_o 

l  .  .  * 

a  i .  .  +  5 .  .  -  is .  . 0 

ij  ij  ij 

(21a) 

,  i 

•  .  i  -i 

a  o  -*•  a 

(21b) 

‘jo 

J 

r1 .  * 

-  i 

-t 

(21c) 

OJ 

J 

rik- 

0 

( 2  Id ) 

r°  ■ 

r°.  >  o 

(21e) 

oo 

01 

4 

r:  - 
10 

r1  *  o 

oo 

(21f  ) 

6 


v'  *.  %  1%  mV  .V 


T 


,1 


MOW 


The  coupling  of  the  Einstein  tensor  in  the  Einstein-Cartan  theory  with 
the  Ray-Smalley  Improved  SEMT  involves  writing  the  field  equations  in  the  so- 
called  self-consistent  form; 


5aB  .  ’  «6u  Suo  a8 

4  Ray-Smalley 


where 


7  (  }  “  V  )  +  2S  a(  ) 

v  \)  \)ct 


(23a) 


and  Q  „  is  the  modified  torsion,  defined  as 
aB 

q  u  =  s  Vi  +  '>6US  v  (23b) 

WctS  aB  "  [a  3]v 

Because  the  mass-conserving  case  is  being  considered,  the  form  of  the  torsion 
in  this  model  becomes 


Q  u  *  5  u  =  S  A 

aB  id  aB 


(23c) 


and  V  (  )  =  7  (  ). 
v  v 


The  field  equations  generated  in  the  model  versus  the  values  of  a  and  3, 
are  summarized  in  Table  3. 
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*_>1 


TABLE  3.  Field  Equations 


GaS  -  7*  (Qa0U  -  Q@ua  • 
a  3  Field  Equations 


Ray-Smalley 


a‘  -  (  jo .  .  a .  .  =  <p(  1  +  e  ) 
ij  iJ 

-  0  i  -  1,2,3 

*  0  i  ■  1,2,3 

2a  -  3(a)2  -  ia„  5.^  +  +  [a,  x]n  »  cP  +  <px21 

12  +  3  ®12  a  +  ^ ^  12  "  ^12°  +  3  a  Si2°^  *  0 

12  +  3  °13  *  +  [3,  tl3  -  iS12°  t23  -  i<px23s12 

21  +  3  521  4  +  [5,  t]21  -  J(S21°  +  3  4  S21°]  *  0 

2  a  -  3(a)2  -  i 5  oij  ■+•  +•  [5,  x]_0  *  <?  +  <px 

i  J  *•  ••  4mm, 

'  ,  r  “  -  l  1^.0-  i  -  21 


21  ?1  L  ’  J21  21  21  J  a  S21  J 
2  a  -  3(a)2  -  ia  oij  ■+•  oin  +  [5,  x],0  *  <?  + 

23  +  3  o23  a  +  [5,  x]23  -  |S21°  ?13  »  iicpx^s21 

i,  +  3  a,,  a  +  [o,  x]_,  -  £S°  t_ _  *  WpT„„s1- 


The  equation  [00],  at  first  appearance,  seems  to  contain  no  spin  energy 
terms.  However,  the  energy  term  is  there  by  virtue  of  the  term,  e,  located  on 
the  right  hand  side  of  the  equation.  Recalling  the  thermodynamic  laws  of  the 
fluid,  as  presented  by  Ray  and  Smalley  [7],  the  differential  of  energy  de  is 
given  by: 


dc  ■  T  ds  *  P  d(l/p)  +  Jw  dsc 

aS 


Although  the  specific  form  of  each  of  the  thermodynamic  variables  for  the 
fluid  is  unknown  in  this  model,  it  may  be  stated  that  the  integrated  energy 
e  represents  a  correction  to  the  standard  energy  term  usually  written  in  the 
standard  theory  of  perfect  fluids.  The  usual  term  for  the  T00  energy  density 
components  of  the  SEHT  is  simply  <p.  In  this  model,  it  is  assumed  that  the 
correction  term  e  is  small  compared  to  the  total  energy  due  to  fluid  density. 
Thus,  in  the  scaling  laws  developed  in  the  next  section,  p*  ■  p  (1  +  e)  will 
have  the  same  scaling  behavior  as  p. 


III.  BEHAVIOR  OF  THE  MODEL  WITH  TORSION 

From  the  field  equations  contained  in  Table  3,  the  following  observa¬ 
tions  are  made.  The  components  [00],  [Oi],  and  [iO]  are  identical  to  their  GR 
counterparts.  Torsion  appears  in  all  the  other  field  equations.  In  order  to 
determine  the  effect  of  torsion  on  the  solutions  to  the  field  equations,  self- 
consistency  must  be  demanded.  To  elucidate,  consider  equations  [12]  and  [21]. 
These  two  are  related  since  the  GR  terms  involve  shear  terms  which  are  sym¬ 
metric,  the  expansion  which  is  a  scalar,  and  the  antisymmetric  product  of 
shear  and  antisymmetric  product  of  shear  and  vorticity.  By  examining  and 
writing  out  equations  [12]  and  [21]  in  more  detail,  they  become 

[12]:  512  “  3  512  *  +  [5>  i]12 

‘  HS12°  +  3isi2°]  =  0  (25a) 

[21]:  a21  +  3  a21  a  +  [5,  t]?1 

-  J[S2l°  +  3  i  S21°]  =■  0  (25b) 

Reversing  the  indices  in  Eq.  (25b),  an  equation  identical  to  Eq-  (25a) 
is  obtained, ^except  for  the  terms  involving  torsion.  (It  can  be  shown  that 
[a,  t ] ^ j  *  [a,  r]^.)  Doing  this  operation,  Eq .  (25b)  becomes 

In)i  S12  +  3  512  *  +  [5>  ;i12 

+  itsi2°  +  3  a  Su°]  *  0  (26) 

Comparing  Eq.  (26)  with  (25a),  it  is  immediately  seen  that  the  term 
involving  torsion  (fourth  term  in  brackets)  must  be  zero.  Thus,  the  fourth 
term  must  satisfy  the  identity 

[S12°  +  3  a  S  12° ]  -  0  (27) 

which  allows  us  to  solve  for  the  torsion  as  a  function  of  time: 

S12°  =  S12°(0)  e'3a  •  (28) 

Examining  field  equation  [13],  an  identity  can  immediately  be  made  by 
recalling  that 
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O  ,  12  o  .  12 

u  »  s  u  *  2<ps 


(29) 


m 


% 


pve 

vs 


r>,'  •: 

V?.- 

■  V-, 


kj; 

ft! 


i<p 


s 
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since  u°  -  1.  Thus,  substituting  Eq.  (29)  into  the  right-hand  side  of  field 
equation  [13]  and  subtracting  from  both  sides,  it  can  be  written  as 


ai3  +  3  5L2  a  +  [5>  ?1i3  *  (3/2^  s12°  *23  =  0 
while  equation  [31],  using  the  same  identity,  can  be  written 


31 


+  3  53i  “  +  t5*  ^21  ‘  *S  °  ' 


12  k23 


=  0 


(30a) 


(30b) 


Since  the  shear  and  shear-vorticity  commutator  terms  are  symmetric  in  their 
indices,  Eq.  (30b)  can  be  transformed  into  the  following,  by  switching 
indices : 


al3  +  3  a13  a  +  ta*  i]i3  '  *s12°  *23  *  0  (30c) 

To  have  consistency  between  Eqs.  (30a)  and  (30c),  set  t23  *  0.  Similar 
reasoning  using  field  equations  [23]  and  [32]  yields  t(3  *  0. 

Adding  equations^ 11] ,  [22],  and  [33]  together  and  using  the  fact  that 
Trace  [a]  ■  0  and  [a,  t]h  *  0,  the  resulting  equation  can  be  written  as 

*■} 

-6  i  -  9  (a)-  -  (3/2)  5:  *  3<?  +  2 .<o  t  s~~  (31) 

i  . 

Field  equation  [00]  is  written  as 

3  (a)-  -  io  a13  *  <p  (1  +  e),  (32) 

where,  recalling  from  previous  argument  [p1  =p(l+e)],  p  (1  +  e)  can  simply 
be  replaced  by  p  in  the  following  equations.  Thus,  the  time  rate  of  change  of 
the  expansion  can  be  solved  as 

(i)2  -  (1/6)  S:  +  ( 1/3 )tcp  .  (33) 

Equation  (33)  can  be  substituted  in  Eq .  (31)  to  yield 

a  -  -ko1  -  i<(p  +  p)  +  (</3 )  T  s12  (34) 


The  expansion  term  a  can  be  written  in  terras  of  an  apparent  rate  of 
change  of  radius  of  the  universe  (Hubble  expansion)  as 

a  «  R/R 


a  *  R/R  -  (R/R)- 

Substituting  the  above  plus  writing  the  shear  scalar  as  p0 
(34)  can  be  written  as  (let  <  -  1); 

10 


(1/2)  a2,  Eq. 


(35a) 

(35b) 


t 


R/R  *  -  (2/3)  3  -  3/6  -  P/2  +•  (.1/3)  3 


To  arrive  at  the  appropriate  power  laws  for  each  of  the  terms  In 
Eq.  (36),  It  Is  necessary  to  make  use  of  several  Identities.  Using  the  con¬ 
tracted  Blanchi  Identities  [9),  It  Is  possible  to  arrive  at  the  first-order 
differential  equation  relating  density  and  pressure  (letting  P  -  py).  Then 


3  *  -3  a  (3  +  P) 

*  -3  a  ( 1  +  y )  p 


(37a) 
r  37b ) 


gives, 


3  *  3  e 

o 


-  3a(  1-cy  ) 


From  Eq.  (35a),  exp(a)  *  R  can  be  written  and  then  substituted  In  Eq .  (38)  to 
find 


0*3  R 
o 


*3( 1+Y  ) 


Because  of  the  relationship  between  p  and  P,  a  power  law  can  Immediately  be 
written  for  pressure  P  as 


P  -  b  30  R 


-3(1+y) 


The  scaling  for  the  shear  scalar  term,  pa  *  l/ZoijO1^,  can  be  determined  bv 
the  following  analysis.  From  the  form  of  the  torsion  chosen  (spin  vector 
aligned  along  the  z-axls),  the  non-zero  components  of  the  coordinate  shear  Is 
expressed  in  terms  of  a  "shear  ellipsoid"  with  the  tetrads  aligned  such  that 

[9] 


5, , 

( 

1  *  coso)  + 

( 1 

-  coso) 

i  l 

-  * 

— 

X 

( 

.0 

.  -  COSO  )  -  , 

1 1 

-  coso) 

- 

.0 

°33 

"33 

al2 

* 

"  ^12 

coso  +  (;°  - 

:22} 

sino 

(41a) 
(41b) 
(41c) 
(  4  Id  ) 


where 


Z°.  -  (1/2)  5°.  e'3a 
ij 


( 4  le ) 


3*2/  t12  dt 


(  4 1  f  ) 


The  shear  scalar  can  be  represented  as  (sum  of  non-zero  terms) 


°11  +  522  +  533  *  2  512  *  2(511  *  522  +  512  *  511  522) 


/UY 


Performing  the  Indicated  operations  obtains 


, , , .  _o  ,2  ,  _o  .2  ,_o  .2  ,  ,o  „o 

^o  "  ^  “11  +  “22  +  “12  +  L11  “22 


_  o  ,  „  o  _  o  .  .  ,  , 

♦  -  sm2#  j 

In  Eq .  (43),  replace  all  2®  by  their  corresponding  terms, 

.  -o  :  -o  2  -o  :  -o-o 

J  *  i  ,  l  »  i  ;  _  „  -  1  c  .  _  .1  •  j  .  .  z  ,  „ 


-o  -o  -o  "M  C44'l 

*  j,  -  s;n  Ij.e  ' 

Thus,  the  term  o  scales  as  R  **. 

j 

To  determine  the  scaling  laws  for  the  shear  term  ri2  Involved  in 
Eq .  (36),  use  the  result  of  Reference  9  which  shows  by  dimensional  analysis 
that  the  shear  evolution  is  related  to  the  proper  time  T  by  the  relationship 

;  ,  T-1  (45) 

ij  ij 

To  get  a  scaling  relationship  between  R  and  time  T,  proceed  as  follows: 

The  field  equation  [00]  can  be  written  as 

( R/R ) “  -  (1/3)  3°  R'6  *  o  R*3< 1 "y )  (46) 

a  o 

For  y  ■  0  (dust),  the  last  term  is  R-^,  for  y  m  1/3  (radiation)  it  is  R-14,  anc 
for  y  ■  1  (stiff  matter)  it  is  R-^.  Depending  on  the  region  of  Interest  (R  > 
or  <  0)  and  the  value  of  y,  one  or  the  other  term  in  Eq .  (46)  will  dominate. 

As  an  example,  for  dust  ( y  ■  0)  and  R  >  1,  the  second  term  will  dominate. 

Thus,  the  differential  equation  in  Eq .  (46)  is  approximated  and  written  as 

(S/R)2  -  ao  R'3  .  (47) 

Solving  for  R,  the  following  is  obtained: 


)/2»  p 


)-!  •»'  - 


(48a) 


which  indicates  that  R  scales  with  respect  to  time  as 


(48b) 


Proceeding  similarly  for  the  other  values  of  y,  both  for  R  >  1  and  R  <  1,  siml 
lar  scaling  relationships  are  derived  between  R  and  T  which  are  summarized  in 
Table  4. 


% 


TABLE  4.  Scaling  Relationships 


Type  Natter 

Dominant  Term 

R  Values 

Scaling  Law 

1/3  o°  R'6 
a 

R  <  1 

R  %  T1/3 

Dust:  y  *  0 

o°  R*3 

R  >  1 

->  /  -5 

R  ^  T~' 

1/3  If6 

R  <  1 

Radiation:  y 

-  1/3 

o°  R*4 

R  >  1 

R  ^ 

1/3  (0°  +  3  )R'6 

1  0 

R  <  1 

R  %  T*y J 

Stiff  Matter: 

y  ■  1 

1/3  to3  «■  o  R-° 

Z  0 

R  >  1 

R  T*  ^ 

1/3  3 

As  seen,  for  all  matter  types  and  R  <  1,  R  ^  T  .  Thus,  T  goes  as  R  and  the 
scaling  law  for  shear  is  written  as 


To  complete  the  scaling  laws  for  all  the  terms  involved  In  Eq .  (36),  the 
scaling  properties  for  the  spin  density  s^  need  to  be  postulated.  In  order 
to  do  this,  the  first-order  differential  equation  relating  the  time  evolution 
of  spin  density  to  tetrad  rotation  rate  and  as  spin  density  Is  written  [7] 

D( s  .  . )  »  w . .  sl .  +  w . .  sl ,  (50) 

ij  ii  j  jl  i 

which  shows  that  the  time  derivative  of  spin  density  is  proportional  to  terms 
like  w n  s*j.  Thus,  a  differential  equation  relating  the  time  rate  of  spin 
and  spin  and  tetrad  rotation  may  be  written  as 


s 


V  s  , 


(51a) 


1  3 


which  gives  (w 


-T 


t*/T) 


ds/s  ■  -t°  dT/T  (51b) 

Solving  for  s,  gives 

Sp  ^  1/T  ^  R'3  (51c) 

Combining  the  scaling  laws  for  density,  shear,  and  spin,  the  following  is 
obtained : 

(1/3)  p  Tp  S12  -  (1/3)  oo  T°;  si2R*3(3^}  ( 5  Id ) 

Combining  Eqs.  (39),  (40),  (44),  and  ( 5 Id )  In  Eq .  (36),  derives 

R/R  -  -(y/2  +  1/b)  o  R'3(i^}  -  C/3)p°  R'6 

o  3 


+ 


(1/3 )  P 

o 


r-3(3^) 


(52) 


The  usual  relationships  between  matter  density  and  pressure  used  In  cosmologi¬ 
cal  models  are  non-lnteractlng  dust  (y  ■  0),  radiation  (y  -  1/3),  and  so- 
called  "stiff"  matter  (y  ■  1).  For  dust,  radiation,  and  stiff  matter,  the 
term  containing  shear  and  spin  density  scales  as  R_  ,  R-1®,  and  R”^,  respec¬ 
tively,  Is  Immediately  concluded.  By  putting  in  the  respective  values  for 
y  Into  Eq.  (52a),  the  following  Is  obtained  for  dust  (y  -  0): 

R  -  -(1/3)  pq  R*2  -  C/3)  o°  R*5  *  (1/3)  pq  t°:  s°;  R-9  (53a) 

for  radiation  (y  *  1/3): 


R  -  -(1/3)  p  R 
o 


C/3) 


(  1  /  3  ) 


and,  for  stiff  matter  (y  •  1): 


H  -  -C/3)  (p  -  o°)R° 

O  3 


C  /  3  )  o 

o 


(53b) 


(53c) 


As  seen  from  Eqs.  (53a,  b,  and  c),  for  small  values  of  R  (less  than  one  In 
normalized  coordinates),  the  spin  density  term  will  dominate  as  a  large  posi¬ 
tive  term,  thus  providing  the  source  for  Increasing  expansion;  shear  does  not 
effectively  dampen  the  expansion.  According  to  this  model,  a  radiation  domi¬ 
nated,  early  epoch  cosmology  will  have  its  expansion  strongly  driven  bv  the 
spin  kinetic  energy  of  the  fluid. 


IV.  CONCLUSIONS 


By  exact  solution.  It  Is  shown  that  expansion.  In  the  early  radiation 
doalnated  phases  of  a  Blanchl  Type  I  Elnsteln-Cartan  cosaology,  is  driven 
positively  by  the  spin  kinetic  energy  of  a  perfect  fluid  using  the  improved 
SENT  of  Ray  and  Saalley.  It  is  immediately  obvious  that  R  ■  0  is  not  a  solu¬ 
tion  to  either  Eq.  (53a,  b,  or  c).  A  solution  set  to  these  equations  can  be 
found,  but  their  exact  fora  depends  on  the  boundary  conditions.  Even  in  this 
slaple  aodel,  It  is  concluded  that  the  spin-energy  not  only  drives  the  expan¬ 
sion,  but  prevents  the  occurrence  of  the  singularity  and  causes  the  apparent 
radius  of  the  universe  to  “bounce". 

Gasperlni  [5]  has  deaonstrated  a  spin-driven  inflation  using  a  time- 
averaging  and  scaling  analysis  of  the  Elnsteln-Cartan  equations.  Kopczynskl 
[10]  and  Trautaan  [11]  have  obtained  minimal  radius  solutions  for  torsion 
cosmological  models  containing  polarized  dust. 

When  going  to  an  Elnsteln-Cartan  cosmology  with  spin  density  using  the 
Ray-Smalley  Improved  SEMT  of  spinning  fluids,  It  Is  apparent  from  this  simple 
model  that  the  properties  of  the  standard  Blanchl  Type  I  cosmology  are  dras¬ 
tically  changed.  Thus,  it  could  be  argued  that  astronomical  observations 
which  lead  one  to  classify  behavior  as  a  Blanchl  Type  of  higher  number  In 
standard  general  relativistic  theory,  could,  in  reality,  be  torsion  "masquer¬ 
ading"  Itself  as  some  sort  of  pseudo-curvature  in  a  universe  which  obeys  the 
Elnsteln-Cartan  formalism.  Much  work  remains  to  be  done  in  this  area,  in¬ 
cluding  a  systematic  re-examination  of  cosmological  properties  of  all  Blanchl 
Type  structures  within  the  framework  of  the  correct  description  of  spinning 
fluids  In  the  torsion  theory  of  gravitation. 
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